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Most  of  the  past  work  on  axi symmetric  wake  flows  has  dealt 

with  simple  shapes  such  as  circular  disks,  ellipsoids,  spheroids 

and  some  slender  bodies.  The  vast  majority  of  these  studies  moreover 

has  dealt  with  unpropelled  configurations. 

The  work  on  propelled  bodies  has  been  primarily  restricted  to 

jets  injected  from  circular  disks  such  as  the  works  by  Ridjanovic1 
2 

and  Haudascher  . Only  in  the  last  four  years  has  the  wake  of 
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propel ler-driven  bodies  been  considered.  Gran  in  1972  and  1973 
studied  the  wake  of  a propeller-driven  Rankine  Ovoid  at  Reynolds 

4 

numbers  based  on  diameter  of  « 6 x 10  . More  extensive  work  on 

4 

propeller-driven  slender  bodies  was  published  in  1974  by  Swanson  , 

5 5 

et.  al . , and  Chieng  , et.  al.,  at  Re^  of  = 6.18  x 10  . From  these 

studies  it  was  found  that  dramatic  changes  occur  in  the  wakes  of 

non-propelled  slender  bodies  and  self-propelled  bodies  of  the  same 

configuration.  Table  I gives  a survey  of  past  wake  studies  on  axi- 

symmetric  bodies. 

In  this  report  several  previously  unstudied  effects  on  the  turbulent 
wake  of  propeller-driven  bodies  are  examined.  The  first  is  the  effect 
of  body  pitch  angle  on  the  wake  properties.  The  second  is  the  addition  of 
an  appendage  such  as  a "sail"  to  an  axisymmetric  body.  The  third  is  the 


evray  1968  18  x x 6 x 1 Spher- 
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effect  of  replacing  a single  propeller  by  an  equivalent  set  of  side- 
by-side  counter  rotating  propellers.  Two  supporting  studies  were  also 
conducted.  One  dealt  with  measurements  of  the  aero-hydrodynamic  forces 
and  movements  on  propeller-driven  slender  bodies.  The  second  involved 
boundary  layer  measurements  on  the  body  itself. 

For  this  work,  a realistic  body  shape  was  chosen.  A model 
fineness  ratio  of  12:1  was  used  because  of  its  similarity  to  most  modern 
submarines,  and  the  appendage  was  scaled  to  simulate  a submarine  sail. 

To  study  the  effect  of  pitch  angle  on  the  wake  properties  the  tests  were 
conducted  at  model  pitch  angles  of  0°  and  -2°.  The  dual -propel  1 er  model 
had  the  same  fore-  and  middle-body  shape  and  the  same  appendage. 

The  primary  testing  was  conducted  in  the  VPI  six-foot,  subsonic 

5 

wind  tunnel  at  Re^  - 4.4  x 10  . This  wind  tunnel  has  proven  in  the  past 
to  be  an  excellent  facility  for  detailed  flowfield  and  turbulence  measure- 
ments because  of  its  low  free-stream  turbulence  level  and  its  uniform 
pressure  field.  Measurements  of  total  pressure,  static  pressure,  flow 
angularity  and  axial  turbulence  intensity  were  made  at  discrete  points 
across  the  wake  at  three  downstream  stations;  2,  10,  and  40  model  diameters 
and  on  the  body  at  several  stations.  These  measurements  then  yielded  all 
velocity  components  and  axial  turbulence  intensity  through  the  wake. 


fh 

The  measurements  were  taken  using  a yawhead  pressure  probe, l?pitot- 
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static  probe,  a straight  hot-wire  sensor  and  a boundary  lay(|*  pitot 
rake. 

This  report  is  arranged  into  several  major  sections.  H;rst, 
the  apparatus,  models  and  experimental  methods  are  described  fR  detail 
The  results  are  presented  in  separate  sections,  starting  with  the  force 
and  moment  data.  The  next  subsections  in  the  Resul ts  section  deal 
with  the  body  boundary  layer  measurements,  the  single  propeller  model 
wake  results  and  finally  the  double-propeller  model  wake  results. 

The  results  are  presented  mainly  in  graphical  form,  but  the  data 
for  the  main  test  program  with  the  single  propeller  model  is  also 
tabulated.  Comparisons  with  previous  work  on  propeller-driven  bodies 
are  given  whenever  available. 
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II.  EXPERIMENTAL  APPARATUS  AND  INSTRUMENTATION 


Facility 

All  tests  were  conducted  in  the  VPI&SU  6'  x 6'  subsonic  stability 
tunnel  at  a baseline  condition  of  a dynamic  pressure  of  5.0  inches  of 
water  (approximately  157  ft. /sec.)  yielding  a Re^  based  on  diameter 

5 

4.4  x 10  . The  facility  is  a continuous,  closed  jet,  single  return 
wind  tunnel.  Its  test  section  is  28  ft.  long  allowing  wake  measurements 
at  40  model  diameters  with  the  Model  located  at  the  far  upstream  end  of 
the  test  section.  The  air  stream  has  a low  turbulence  factor  of  1.08. 

To  control  the  tunnel  velocity,  a Barocel  Electric  Manometer  is  used  to 
read  dynamic  pressure  from  a pitot-static  tube  mounted  on  the  test  section 
wall  out  of  the  wake  of  the  model.  Free-stream  temperature  is  monitored 
by  a temperature  probe  on  the  test  section  wall  and  recorded  on  a Digitemp 
temperature  gauge.  Static  free-stream  pressure  is  measured  by  a Validyne 
digital  barometer  Model  DB99. 

Models 

The  model  considered  in  this  investigation  has  an  overall  length 
of  72  inches  and  a maximum  diameter  of  6 inches,  giving  a fineness  ratio 
of  12:1.  The  model  was  strut  mounted  from  the  ceiling  of  the  wind  tunnel 
just  forward  of  the  test  section  (Fig.  1).  The  forebody  is  parabolic  and 
machined  from  laminated  layers  of  plexiglas.  The  centerbody  is  an  aluminum 
cylindrical  tube  with  a simulated  sail,  and  the  tail  bodies  are  also  plexi- 
glas. For  the  overall  dimensions  of  the  body,  see  Fig.  2. 


Both  tail  bodies  were  ogives.  The  first  body  considered  had  a 
single  shaft  which  extended  out  the  stern.  The  second  body  had  twin  shafts 
which  flared  out  at  a 7°  angle  through  the  stern.  Both  models  used  2.75 
h.p.,  DC  motors  to  drive  the  propeller  shafts.  The  single  shaft  model  was 
direct  drive,  while  the  dual  shaft  model  used  a belt  driven  1:1.5  step-up 
from  the  motor  to  the  counterrotating  shafts  (Fig.  3a).  This  was  necessary 
to  achieve  the  higher  rpm's  needed  for  this  model  in  the  self-propelled  mode 
The  single  shaft  model  used  a 6"  diameter,  3 bladed  model  airplane 
propeller,  while  the  dual  shaft  model  used  two  4.375"  diameter,  3 bladed 
propellers  which  were  cut  down  versions  (RH  and  LH)  of  the  single  pro- 
pellers. All  propellers  were  heated  at  the  root  and  twisted  to  a higher 
pitch  to  operate  more  efficiently  at  the  high  air  speeds  used.  The 
effective  pitch  of  the  larger  propellers  as  modified  was  2.46  (Ref.  (4)). 

The  smaller  propellers  were  twisted  to  the  same  angle. 


• *t 


Force  and  Moment  Instrumentation 

The  forces  and  moments  were  measured  by  a six  component  strain  gage 
balance  made  by  Modern  Machine  Tool  Company  internally  mounted  at  the 
c.  g.  of  the  models.  The  arrangement  is  shown  in  Fig.  4.  Note  that  the 
mounting  strut  is  a thin  (0.5  in.)  blade  that  clears  the  "sail"  on  the 
model.  Forces  on  the  sail,  but  not  on  the  strut,  are  measured  along  with 
those  on  the  main  body.  Output  from  the  balance  was  read  on  a Doric 
digital  voltmeter.  Other  support  equipment  included  a Beckman  Counter 
Model  2D-3  which  gave  continuous  readings  from  a magnetic  pickup 
attached  near  the  propeller  shafts  to  determine  the  rpm  of  the  propellers 
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for  the  self-propelled  tests.  To  power  the  d.c.  motors  in  both  models, 
a Sorenson  150  volt,  15  amp  power  supply  was  used.  Pressure  readings  to 
determine  tunnel  speed  were  made  with  a Validyne  digital  barometer  Model 
DB99. 

Make  Pressure  Probes 

The  mean  flow  measurements  were  made  using  a three-dimensional 
Yawhead  probe  constructed  by  United  Sensor.  The  probe  was  used  in  the 
wake  region  when  substantial  flow  angularity  occurred  (Fig.  5).  In  the 
far  wake,  the  flow  angularity  was  sufficiently  small  to  use  a standard 
straight  pitot-static  tube. 

The  output  from  the  pressure  probes  was  measured  on  a Barocel 
Electronic  Manometer  and  displayed  on  a Doric  Digital  Voltmeter.  Graph- 
ical displays  were  provided  by  a HP  7100B  Strip  Chart  Recorder  and  a 
Model  20-3  X-Y  Recorder.  To  read  all  five  pressure  ports  from  the 
yawhead  probe  on  the  electronic  manometer,  a Scanivalve  Model  W0601 / IP-RT 
fluid  switch  wafer  with  solenoid  drive  was  used. 

Hot  Wire  Instrumentation 

The  axial  turbulence  data  was  obtained  using  a straight  hot  wire 
(TSI  Model  1210).  This  sensor  uses  a platinum  plated  tungsten  wire 
.00015  inches  in  diameter.  The  probe  was  operated  at  an  overheat  ratio 
of  1.8  using  a constant  temperature  anemometer  module  (TSI  Model  1050) 
and  a Power  Supply  and  Monitor  (TSI  Model  1051-6).  Readings  were  taken 
from  a DISA  55D35  RMS  Meter. 
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Traverse  Mount 


Both  pressure  and  hot  wire  probes  were  moved  through  the  wake 
using  the  traverse  shown  in  Fig.  6.  This  mount  had  a vertical  movement 
of  4 ft.  and  a horizontal  range  of  11.0  in.  It  was  therefore  necessary 
in  the  far  wake  studies  to  move  the  mount  to  obtain  a wider  horizontal 
range.  Both  vertical  and  horizontal  movement  was  accomplished  with 
variable  speed  motors  which  were  controlled  outside  the  test  section. 

The  probe  position  was  monitored  using  a series  of  potentiometers. 

Boundary  Layer  Rake 

The  body  boundary  layer  measurements  were  made  with  a rake  of  pitot 
tubes  as  shown  in  Fig.  7.  Closer  spacing  of  the  tubes  nearest  the  wall 
was  employed  to  produce  better  resolution  of  the  turbulent  boundary 
layer  profiles  anticipated. 

Since  the  rake  was  built  to  allow  mounting  by  straps  around  the  body 
at  virtually  any  axial  station  and  at  variable  peripheral  locations,  a 
pitot-static  probe  was  attached  to  the  outside  of  the  rake.  This  was  to 
permit  static  pressure  determination  wherever  the  rake  was  located.  This 
caused  concern  as  to  possible  interference  effects.  Previous  work  in 
this  area  by  Krause,  Ref.  (22),  indicates  that  the  change  in  static  pressure 
reading  obtained  from  a pitot-static  tube  due  to  an  adjacently  positioned 
tube,  can  be  held  to  a minimum  if  the  spacing  is  approximately  5 diameters. 
In  cases  where  large  gradients  necessitate  placing  total -pressure  tubes 
closer  than  5 diameters  to  static-pressure  tubes,  the  proximity  effects 
can  be  minimized  by  proper  orientation  of  the  leading  edge  of  adjacent 
tubes  in  relation  to  the  static  orifices.  A logical  extension  of  these 
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results  was  felt  to  be  that  adjacent  total -pressure  tubes,  sans 
static  ports,  should  also  be  aligned  at  their  leading  edge.  Still,  the 
proximity  effects  do  increase  with  closer  spacing  and  thus  a minimum 
tube  spacing  of  2.5  diameters  was  specified  for  the  purposes  of  this 
investigation. 

The  literature  indicates,  also,  that  the  ratio  of  the  distance 

between  the  pitot-static  tube  leading  edge  and  the  static  ports  to  the 

23 

distance  between  those  ports  and  the  support  structure  is  critical 
Depending  upon  whether  the  design  is  "standard"  or  according  to  Prandtl , 
the  recommended  ratio  is  either  1:2  or  1:3.  In  the  standard  design 
the  tube  is  bent  90°  in  a curved  fashion,  16  diameters  behind  the 
static  ports,  whereas  in  the  Prandtl  design  this  bend  is  squared  off 
8-10  diameters  behind  the  ports. 

Due  to  the  nature  of  the  rake  design,  a support  strut  extended  from 
the  base  up  past  the  highest  total -pressure  tube.  By  affixing  a 1/1 6th 
inch  diameter  pitot-static  tube  to  the  top  of  this,  a Prandtl  config- 
uration was,  in  effect,  achieved  with  the  strut  serving  the  function  of  a 
squared  off  tube  bend.  Therefore,  the  pitot-static  tube  was  positioned 
such  that  its  static  orifices  lay  10  diameters  forward  of  the  strut. 
Unfortunately,  we  were  then  faced  with  the  fact  that  the  tip  ends  of  the 
underlying  total -pressure  tubes  were  not  in  line  with  the  pitot-static 
tubes  static  orifices,  as  per  Krauses'  recommendation  (Ref.  (22))  and 
thus  a test  of  this  new  configuration  was  performed  in  the  VP  I open- 
throat,  return-type  wind  tunnel  with  a three  foot  diameter  test  section. 

A four-foot  length  of  2-1/4  inches  diameter  plastic  pipe  was  rigidly 
suspended  by  wire  in  the  test  section  and  fitted  with  a smooth  nose-cone. 
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The  rake  was  attached  to  this  at  various  downstream  positions  and 
connected  to  all  the  same  electronic  equipment  as  was  used  during  test 
runs  on  the  submarine  model  in  the  VPI  Stability  wind  tunnel.  The 
model  was  tested  at.  several  velocities  up  to  3.5  inches  of  water,  and 
static  pressure  data  were  taken  at  selected  points  on  the  model  with 

i 

the  rake.  At  these  same  points,  static  pressure  data  were  obtained 
using  a solitary  pitot-static  tube  suspended  to  the  same  spatial 
location  as  the  rake,  pitot-static  tube  previously  was.  From  these 
data,  plots  of  error  in  static  reading  as  a function  of  both  velocity, 
boundary  layer  height,  and  xVd»  the  nondimensional ized  variable  denoting 
distance  back  from  the  leading  edge,  were  made.  It  can  be  seen  in  Figure 
8 that  the  data  for  2.5  and  3.0  inches  of  water  overlay  each  other  quite 
well  and  indicate  that  the  error  ranges  from  near  5%,  when  the  boundary 
layer  is  very  small,  to  about  21.  The  assumption  was  made  that  the 
curve,  6 versus  l error,  did  not  vary  with  velocity.  This  correction  was 
applied  to  the  submarine  model  measurements. 


III.  EXPERIMENTAL  METHODS 


Force  and  Moment  Tests  and  Self-Propelled  Status 

With  the  model  mounted  in  the  tunnel,  potentiometers  were  used  to 
zero  the  output  from  the  balance  due  to  the  weight  of  the  model. 

The  free-stream  velocity  was  increased  to  5.0  inches  of  water 
while  the  axial  force  was  monitored.  The  propeller  rpm  was  then 
increased  until  the  axial  force  reading  was  zero.  This  was  done  for 
several  angles  of  attack  between  -6°  and  +6°.  At  each  angle  the  pro- 
peller RPM  needed  to  zero  the  axial  force  was  recorded. 

To  expedite  data  reduction,  the  calibration  curves  for  the  six  com- 
ponent balance  were  programmed  in  a data  reduction  routine  for  use  on  the 
VPI  IBM/370  computer.  The  program  converted  the  strain  gage  output  in 
millivolts  into  body  forces  and  moments  in  standard  English  units.  These 
forces  were  then  non-dimensionalized  into  the  standard  body  force  and 
moment  coefficients  using  the  free-stream  static  pressure  and  tempera- 
ture. The  program  also  performs  the  axes  transformation  to  obtain 

and  CD  in  the  wind  axes.  The  results  were  then  machine  plotted  using  an 
interactive  data  reduction  and  plotting  routine. 

Mean  Flow  in  the  Wake 

The  mean  flow  velocities  as  well  as  flow  angularity  were  determined 
using  a five-port  yawhead  probe  in  regions  of  large  flow  angle.  In  the  far 
wake,  a conventional  pitot-static  tube  was  used  and  data  reduction  was  by 
standard  methods.  The  yawhead  probe  was  calibrated  over  a range  of  flow 
pitch  and  flow  yaw  angles  between  -30°  and  +30°. 

It  was  decided  that  an  adequate  description  of  the  wake  could  be 
achieved  with  three  transverse  cuts.  These  cuts  are  seen  in  Fig.  9.  These 


three  traverses  were  made  at  downstream  stations  of  Z/D  = 2,  10,  and  40 
with  the  model  at  aQ  = 0°  and  ag  = -2°  for  the  single  propeller  model  and 
at  only  ag  = 0 for  the  dual  propeller  model.  Measurements  were  taken  at 
0.5  inch  intervals  across  the  wake. 

To  obtain  the  mean  velocity  components  from  the  pressure  measurements 
the  following  relation  for  total  pressure  was  used: 

p0  = Ps  + I pV  + \ ^ 

where  V is  the  mean  total  velocity  and  V'  is  the  sum  of  the  mean  fluctua- 

2 _ 2 1-2 

tions.  Since  V »V'  ,the  j pV'  was  neglected.  The  yaw-head  probe  gave 
direct  readings  of  - Pstoo>  where  is  the  mean  total  pressure  and 
P$too  is  the  free  stream  static  pressure. 


Therefore  since 


P1  ■ Pst  + ? 2 


the  mean  total  velocity  can  be  calculated  from 


total 


J2«p  1 - pst->  - (pst  - pstJ) 

» p 


I . 


then  the  mean  axial  velocity  can  be  obtained  by 

(Jz  = Vcos(FP)cos(FY) 

where  FP  is  the  flow  pitch  angle  and  FY  is  the  flow  yaw  angle.  Similarly 

the  other  velocity  components  are  given  by 

Ux  = Vcos(FP)sin(FY) 

U = V sin(FP) 

y 

To  obtain  the  local  static  pressure  from  the  yawhead  probe  a method 
described  by  Winternitz  (Ref.  6)  was  used.  Since  this  method  is  only 
good  for  substantial  flow  angularity  in  one  direction,  calibration  curves 
were  obtained  for  both  flow  pitch  and  yaw.  The  flow  yaw  angle  FY  is 
found  from  the  calibration  probe  constant 


. P3  ' P2 
6 " p:  - P4 

for  each  station.  Defining  the  hole  constant  as 

P.  - P 

K,(6)  - 

where  P.  is  the  hole  pressure  and  q is  the  dynamic  pressure.  Wintemitz 
showed  that 

. <P1  • Pst->  • <P2  - Pst-> 


K1  ' K2 


<Pst  - Pst->  ■ <P2  - Pst->  - “V 


Therefore 


<pst  - pst->  " 


K1<P2  - PstJ  - ¥P1  - PstJ 

Kl  - K2 


Similarly  for  the  pitched  case  where  the  yaw  angles  are  small 

K - P4  ~ Ps 
“ ' P1  - P2 


(Pst  ‘ PstJ 


K1(P5  - PstJ  - K5^P1  - PstJ 
K1  - K5 


From  these  relations  the  local  static  pressure  could  be  determined  when 
there  was  substantial  pitch  angles  with  small  yaw  angles  or  for  large  yaw 
angles  with  small  variations  in  pitch. 


Turbulence  Measurements 

The  axial  turbulence  intensity  was  measured  at  the  same  points  as 
the  mean  flow  measurements  using  a straight  hot-wire.  The  constant 
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temperature  method  was  used  as  discussed  in  Ref.  8.  The  hot  wires  were 
calibrated  by  measuring  their  total  voltage  at  different  dynamic  pressures 
between  0 and  8.0  inches  of  water.  This  calibration  was  done  at  three 
free-stream  temperatures,  72°,  84°,  and  92°  F to  see  how  the  calibrations 

were  affected  by  temperature.  By  measuring  the  total  voltage  E and  the 

_ 2 

RMS  voltage  e‘  the  axial  turbulence  could  be  obtained  from  the  relation 

u'2  _ 2E  e* 2 
U ' n(E2  - E02)" 

where  n = .45  for  Re  < 44  based  on  the  diameter  of  the  hot  wire,  and  EQ 
is  the  total  voltage  reading  of  zero  velocity.  Also  by  using  the  cali- 
bration curves  the  mean  axial  velocity  could  be  found  and  compared  with 
the  mean  velocity  obtained  from  the  pressure  probe. 

Boundary  Layer  Measurements 

Measurements  were  made  at  X'/D  = 1,  6,  8,  10  an  d 11.17*  and  at 
peripheral  locations  in  the  line  with  the  sail  (e  = 0),  opposite  the 
sail  (e  = 180°)  and  halfway  around  the  body  (o  = 90°).  The  data  were 
reduced  using  the  Bernoulli  equation,  the  local  pitot  pressure  and  the 
static  pressure  determined  using  the  pitot-static  probe  at  the  top  of 
the  rake. 


* X'/D  for  the  body  boundary  layer  measurements  are  measured  from  the 
body  nose.  Wake  measurements  use  X/D  measured  from  the  stern. 
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IV.  RESULTS 


Force  and  Moment  Data 

These  tests  were  conducted  with  the  models  in  two  locations  in  the 
tunnel  test  section.  Tests  were  run  with  the  models  in  the  very  front 
of  the  test  section  to  determine  the  self-propelled  conditions  for  the 
wake  tests.  Tests  were  also  run  with  the  models  in  the  middle  of  the 
test  section,  where  the  flow  is  more  uniform,  in  order  to  develop  basic 
force  and  moment  data  for  this  body  shape.  The  latter  data  are  presented 
in  this  section. 

The  single  shaft  model  without  a propeller  was  strut  mounted  from 
the  ceiling  of  the  tunnel  at  zero  geometric  pitch  angle.  To  align  the 
model  with  the  free-stream  velocity  vector,  a run  was  made,  the  tunnel 
was  stopped,  adjustments  were  made,  and  tests  were  rerun  until  no  side 
force  was  measured.  After  proper  alignment  was  achieved,  the  six  force 
and  moment  readings  were  recorded  from  the  voltmeter.  Next,  the  angle 
of  the  pitch  was  changed,  and  the  six  readings  were  recorded  again.  This 
procedure  was  repeated  until  measurements  had  been  taken  for  pitch  angles 
between  -6°  and  6°. 

The  next  step  was  to  attach  the  propeller  and  repeat  the  measure- 
ments for  the  self-propelled  case.  To  do  this,  the  wind  tunnel  was  run 
at  test  speed  while  the  axial  force  reading  was  monitored.  Then,  by 
slowly  increasing  the  voltage  to  the  DC  motor  in  the  model,  the  rpm  of 
the  propeller  was  increased  until  the  axial  force  reading  was  zero.  Once 
the  axial  force  reading  was  zeroed,  the  other  five  readings  were  taken, 
and  this  procedure  was  repeated  for  all  pitch  angles. 
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To  study  the  effect  of  twin  propellers  on  the  force  and  moments, the 
single  shaft  stern  section  was  replaced  by  the  twin  shaft  tail  section. 


I 


I 


s 

1 


The  tunnel  was  again  run  at  q = 5.0  in.  H^O  (157  ft./sec.),  and  the  side 
force  was  monitored  with  no  props  to  check  for  proper  alignment.  Once 
proper  alignment  was  established,  the  force  and  moment  readings  were 
recorded  for  the  same  range  of  pitch  angles.  The  props  were  then  attached 
to  the  shafts  and  the  same  procedure  was  followed  as  before  for  the  self- 
propelled  configuration. 

All  forces  and  moments  were  put  in  coefficient  form  and  plotted  ver- 
sus the  angle  of  pitch.  For  each  model  ,the  self-propelled  and  non- 
propelled  cases  were  plotted  together  for  comparison.  There  were  no 
significant  side  forces  or  yawing  moments  measured  with  proper  alignment 
of  the  models.  The  results  are  presented  in  Fig.  10-16. 

Consider  first  the  results  for  the  single  propeller  model.  In  Fig. 

10,  the  drag  coefficient  is  shown  versus  pitch  angle  for  both  the  unpro- 
pelled (no  propeller)  and  self-propelled  cases.  A true  self-propelled 
condition  was  achieved  by  simply  zeroing  the  axial  force  during  the  test 

for  -4°  < a < 4°.  This  procedure  produced  small  errors  ,r  ^ (Cp)unpropel led 

'lD  ~ 10  1 

for  larger  pitch  angles.  The  drag  of  the  unpropelled  body  increased  rather 
sharply  for  positive  pitch  angles. 

The  lift  coefficient  data  are  given  in  Fig.  11  where  it  can  be  seen 
that  the  addition  of  running  propellers  increases  the  lift  curve  slope. 

The  change  is  more  pronounced  at  positive  pitch  angles.  The  increase  in 
lift  curve  slope  is  due  to  at  least  two  effects.  The  first  is  the  compo- 
nent of  thrust  normal  to  the  wind  vector  which  gives  an  increase  in 
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proportional  to  the  sine  of  the  pitch  angle.  The  second  effect  is  a nor- 
mal body  force  on  the  propeller  when  the  body  is  pitched  produced  by 
turning  of  the  fluid  passing  through  it. 

Pitching  moment  results  are  presented  in  Fig.  12.  A small  decrease 
in  the  slope  of  the  curve  is  produced  by  the  propeller.  Since  the  thrust 
is  always  along  the  body  axis,  only  the  effects  of  fluid  turning  by  the 
propeller  when  the  body  is  pitched  are  felt  here. 

Figs.  13,  14,  and  15  contain  the  results  obtained  for  the  dual  propel- 
ler model.  The  drag  coefficient  results  in  Fig.  13  are  qualitatively 
similar  to  those  for  the  single  propeller  case.  The  no-propeller  drag  is 
somewhat  higher  .presumably  due  to  the  presence  of  the  exposed  shafts. 

The  effect  of  propulsion  on  the  variation  of  the  lift  coefficient 
with  pitch  angle  (Fig.  14)  is  perhaps  less  but  similar  in  character  to  that 
on  the  single  propeller  model.  On  the  other  hand,  the  spread  for  the 
propelled  and  unpropelled  arrangements  for  pitching  moment  in  this  case 
is  larger  than  that  for  the  single-propeller  model. 

Lastly,  we  show  rolling  moment  results  for  both  models  in  Fig.  16  in 
the  self-propelled  mode.  The  unbalanced  torque  of  the  single  propeller 
model  can  be  clearly  seen,  while  the  counter-rotating,  dual  propellers 
show  the  expected  zero  net  torque. 

Single  Propeller  Model  Wake  Surveys 

Since  the  sail  of  the  model  as  well  as  the  effect  of  pitch  angle 
was  to  induce  asymmetries  in  the  wake,  it  was  necessary  to  probe  the 
entire  wake  structure.  For  this  reason  all  of  the  profiles  presented 
extend  from  one  edge  of  the  wake  to  the  other  or,  in  the  case  of  the 


vertical  traverse,  from  the  edge  of  the  wake  into  the  wake  region  of 
the  sail.  The  mean  axial  velocity,  static  pressure,  flow  angularity  and 
axial  turbulence  intensity  are  plotted  for  both  a&  = 0°  and  -2°  at  stations 
Z/D  = 2,  10,  and  40.  The  coordinate  systems  and  sign  conventions  used  are 
shown  in  Fig.  17. 

First, consider  the  results  for  the  model  at  ag  = 0°.  In  the  near 
wake  (Z/D  = 2)  there  is  a momentum  defect  at  the  center  of  the  wake  and 
then  a momentum  excess  region  in  the  outer  portion  of  the  wake  (Figs.  18 
and  19).  The  static  pressure  decreases  in  the  wake  of  the  propeller  as 
in  the  case  of  a swirling  vortex.  The  propeller  swirl  is  indicated  by 
examining  the  flow  angularity.  For  the  horizontal  profiles  (Fig.  18), 
large  variations  in  flow  pitch  are  measured  with  very  small  changes  in  the 
flow  yaw,  but  in  the  vertical  profiles  (Fig.  19)  the  substantial  flow 
angularity  is  in  the  yaw  direction. 

At  Z/D  = 10  the  momentum  excess  and  defect  regions  have  decreased 
and  spread  out.  The  static  pressure  variation  has  also  decreased  and 
the  swirl  has  decreased  to  the  point  that  the  flow  angularity  is  only 
half  the  value  measured  at  Z/D  = 2 (Fig.  20). 

By  Z/D  = 40  there  is  only  a small  momentum  excess  in  the  wake  and 
its  character  is  very  similar  to  the  far  wake  of  a zero-momentum  jet 
(Refs.  1 and  2).  The  static  pressure  is  uniform  across  the  wake  and 
the  swirl  has  completely  diminished  as  seen  by  the  small  flow  angu- 
larity (Figs.  21  and  22) . 

The  axial  turbulence  intensity  is  characterized  by  a highly 
fluctuating  nature  in  the  near-wake  (Z/D  =2).  At  this  station, the 


velocity  fluctuations  are  up  to  10%  of  the  free-stream  velocity  in 
the  center  of  the  wake  (Fig.  23).  At  Z/D  = 10  the  turbulence  intensity 
has  decreased  to  about  6%  of  the  free-stream  velocity  and  does  not  vary 
significantly  across  the  wake  (Fig.  24).  The  maximum  value  has  fallen 
to  about  2%  by  X/D  = 40  (Fig.  25). 

The  effect  of  the  sail  on  the  mean  flow  and  turbulence  quantities 
is  seen  in  the  vertical  profiles  at  each  station.  The  sail  produces  a 
larger  momentum  deficit  in  the  upper  portion  of  the  vertical  profile 
which  indicates  the  increased  drag  on  the  body  due  to  the  sail  (Figs.  26 
and  27).  The  flow  angularity  is  also  reduced  in  the  region  behind  the 
sail.  This  reduction  is  as  much  as  6.5°  at  Z/D  = 2,  (Fig.  19),  and  3° 
at  Z/D  = 10.  Also  the  static  pressure  is  increased  substantially  by 
the  presence  of  the  sail.  This  increase  is  up  to  50%  in  the  upper 
portion  of  the  vertical  profile  as  compared  to  the  lower  portion  (Fig.  19). 
The  effect  on  axial  turbulence  intensity  was  to  increase  the  level  of 
turbulence  in  the  vertical  profile  behind  the  sail  (Figs.  28  and  29). 

The  mean-flow  results  for  the  pitched  condition  (Figs.  30-33)  show 
increased  asymmetries  in  the  vertical  profiles.  The  flow  angularity 
is  decreased  in  the  upper  region  of  the  vertical  profile  by  9°  at  the 
Z/D  = 2 station  (Fig.  31)  and  by  6°  at  the  Z/D  = 10  station  (Fig.  33). 

The  static  pressure  is  again  higher  above  the  wake  center  for  the  pitched 
configuration.  The  maximum  turbulence  intensity  is  increased  by  15% 
at  the  Z/D  = 2 station  (Fig.  28)  and  by  5%  at  the  Z/D  = 10  station  (Figs. 

34  and  35).  At  40  diameters,  the  turbulence  intensity  and  mean  flow  profiles 
are  about  the  same  for  both  the  pitched  and  unpitched  configurations  (Fig. 

36  and  37).  The  comparison  of  maximum  turbulence  intensity  with  our  previous 
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work  shows  an  increased  turbulence  level  due  to  the  higher  propeller  rpm 
needed  to  balance  the  extra  drag  of  the  sail  (Fig.  38)  compared  to  the 
work  in  Refs.  (4)  and  (5). 

From  the  results  presented, several  conclusions  about  the  effects 
of  appendages  and  non-zero  pitch  angles  on  the  wake  structure  of  a pro- 
peller driven  body  can  be  asserted. 

First , an  appendage  such  as  a submarine  sail  reduces  the  propeller 
swirl  in  the  region  subjected  to  the  wake  of  the  sail.  This  decrease 
in  swirl  causes  an  increase  in  static  pressure  as  would  be  seen  in  a vor- 
tex of  decreasing  strength.  The  axial  turbulence  intensity  is  increased 
by  the  value  of  the  sail  induced  turbulence  at  each  station.  This  sug- 
gests that  the  turbulence  in  the  wake  of  an  appendage  might  be  super- 
imposed on  the  turbulence  profile  of  the  main  body. 

Second  , the  pitched  condition  further  reduces  the  swirl  and  increases 
the  static  pressure  in  the  region  behind  the  sail.  These  effects  are 
primarily  confined  to  the  vertical  plane  and  little  change  in  the 
horizontal  profiles  is  seen.  The  maximum  turbulence  intensity  in  the 
wake  up  to  Z/D  = 10  is  increased  for  the  pitch  condition.  However,  by 
40  diameters  the  turbulence  levels  are  equivalent  for  both  the  pitched 
and  unpitched  conditions. 

Third,  the  effect  of  the  sail  on  the  drag  coefficient  is  to  increase 

it  from  Cp  = 0.09  as  calculated  for  the  body  without  sail  by  Swanson,  et  al 

★ 

(Ref.  (4))  to  C0  = 0.30  with  the  sail.  This  large  increase  in  CD  in- 
dicates the  need  for  streamlining  of  such  appendages  on  slender  axi- 
symmetric  bodies. 

* Both  these  values  correspond  to  the  model  in  the  far  upstream  part  of  the 
test  section  where  some  effects  of  the  contraction  section  persist. 
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Dual  Propeller  Model  Wake  Surveys 

It  should  be  said  clearly  at  the  outset  that  this  test  series  was 
continuously  plagued  by  mechanical  problems  directly  traceable  to  the 
7°  total -angle  flare  between  the  two  counter-rotating  shafts  (see  Fig. 

3).  The  flexible,  constant  rotational  speed  couplings  were  subject  to 
unpredictable  failure  with  a service  life  of  anywhere  between  roughly 
fifteen  minutes  and  two  hours.  We  have  made  repeated  inquiries  to  all 
possible  suppliers  and  have  been  assured  that  the  items  we  are  using 
are  the  best  available  for  the  job.  A coupling  failure  was  also,  unfor- 
tunately, often  accompanied  by  damage  to  one  or  both  of  the  propellers. 

A second  problem  encountered  in  performing  these  tests  was  the 
amazing  sensitivity  of  the  near-wake  profiles  to  small  variations  in  the 
propellers  and/or  the  model.  Even  with  great  care  exercised  in  making 
the  model  and  sets  of  RH/LH  propellers  as  symmetrical  as  possible,  con- 
siderable side-to-side  asymmetries  were  often  encountered.  It  became 
necessary  to  run  systematic  near-wake  profiles  using  various  combinations 
of  propellers  until  reasonable  symmetry  was  achieved.  This  symmetry 
problem  was  actually  so  severe  that  it  has  prompted  us  to  wonder  what  the 
state  of  affairs  is  behind  real  submarines  with  side-by-side  propellers. 
The  symmetry  problem  also  interacted  in  an  unfavorable  way  with  the 
coupling  failure  problem  mentioned  in  the  first  paragraph.  As  we  pro- 
ceeded through  the  test  program,  sets  of  propellers  were  destroyed  and 
had  to  be  replaced.  The  new  sets  were  not  precise  matches  in  terms  of 
near-wake  profiles  to  the  originals.  Thus,  the  test  results  to  be  pre- 
sented below  were  made  with  three  different  sets  of  propellers 
matched  in  terms  of  gross  parameters  such  as  thrust  but  not  in 
scale  detail. 
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Mean  flow  results  are  shown  in  Figs.  39  - 43.  The  horizontal  cut 
through  the  near-wake  (Z/D  = 2)  shown  in  Fig.  39  shows  the  type  of  sym- 
metry that  we  judged  as  acceptable.  Only  with  prior  viewing  of  some  of 
the  unsymmetrical  profiles  that  we  obtained  with  propeller  sets  that 
looked  acceptable  to  the  naked  eye  can  one  objectively  judge  the  results 
in  Fig.  39.  Notice  that,  as  before,  the  propellers  act  mainly  to  produce 
flow  angle  in  only  one  direction. 

A vertical  cut  down  behind  the  sail  and  over  the  tip  of  the  stern 
is  shown  in  Fig.  40.  Here,  the  interaction  between  the  tips  of  the  two 
propellers  produces  a complicated  flow  angle  pattern. 

A vertical  cut  down  over  one  of  the  two  propeller  hubs  is  shown  in 
Fig.  41.  Reasonable  vertical  symmetry  of  profile  is  achieved  here. 

The  mean  flow  results  for  Z/D  = 10  and  40  were  obtained  with  a 
Pitot-Static  probe,  so  that  much  more  detailed  coverage  of  the  flow  in 
a transverse  plane  was  possible.  Those  results  are  shown  in  Figures  42 
and  43  as  profiles  of  U-Ue  overlaid  with  a silhouette  of  the  model  in 
the  background.  These  plots  permit  a comprehensive  picture  of  the  develop- 
ment of  the  total  flow  field.  By  Z/D  = 10,  the  maximum  velocity  variation 
has  fallen  to  approximately  13%  of  the  freestream  speed  from  its  value 
of  approximately  25%  at  Z/D  = 2.  By  Z/D  = 40,  the  maximum  is  further 
reduced  to  6%.  The  corresponding  values  for  the  single  propeller  model  at 

a = 0 were  28%  at  Z/D  = 2,  12%  at  Z/D  = 10,  and  7%  at  Z/D  = 40  indicating 
e 

about  the  same  level  and  rate  of  decrease  for  the  two  models. 

Our  difficulties  with  propeller  failures  and  symmetry  with  this  dual 
propeller  model  forced  only  a limited  effort  with  the  hot-wire  for  this 
case.  A horizontal  and  an  abbreviated  vertical  traverse  for  axial  turbu- 
lence intensity  was  obtained  at  Z/D  = 2.  The  results  are  shown  in  Fig.  44. 


The  observed  maximum  value  of  12%  is  somewhat  higher  than  the  maximum 


for  the  single  propeller  model  with  = 0.  at  the  same  station.  The 
highest  value  was  observed  where  the  tips  of  the  two  side-by-side 
propeller  circles  overlap.  Another  peak  is  observed  at  the  outer  edge 
of  the  propeller  circle  (X/R  = 1.375). 
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Body  Boundary  Layer  Measurements 


The  primary  purpose  of  this  task  was  to  aid  in  interpreting  the  near- 
wake surveys.  Thus,  these  data  were  taken  with  the  body  mounted  in  its 
forward  location  where  the  flow  is  still  adjusting  slightly  from  the  tun- 
nel contraction  section.  The  data  should  then  be  viewed  as  corresponding 
to  the  uniform  flow  over  a body  which  is  slightly  distorted  from  that 
actually  tested.  Nonetheless,  several  interesting  semi-quantitative  and 
qualitative  observations  regarding  the  boundary  layer  development  as  a 
subject  in  itself  can  be  made  on  the  basis  of  the  data  obtained. 

Before  a discussion  of  the  data,  a description  of  the  reference  coor- 
dinate system  is  necessary  for  better  identification  of  the  subsequent 
graphs.  The  origin  is  located  at  the  nose  of  the  submarine  model  with  the 
X'-axis  as  the  axis  of  symmetry  of  the  submarine  (minus  the  sail).  The 
two  coordinates  X ' / D and  0,  locate  the  rake  with  respect  to  the  submarine. 

The  peripheral  angle,  0,  is  measured  from  a value  of  zero  in  the  vertical 
plane  through  the  sail  and  the  body  axis.  Clockwise  looking  upstream  is 
positive.  The  vertical  coordinate,  Y'/D,  is  measured  from  the  local  body 
surface. 

Figures  45  through  47  are  non-dimensional i zed  velocity  profiles  at 
specified  X'/D  locations,  with  three  values  of  0 at  each.  Figures  45  and 
46  for  X'/D  values  of  6 and  8,  respectively,  show  the  effect  of  the  sail 
(0  = 0°)  as  evidenced  by  the  velocity  profile  characteristics  at  0,  90,  and 
180  degrees.  Also,  note  how  closely  the  velocity  profiles  for  90  and  180 
degrees  agree,  which  both  indicates  that  the  flow  disturbed  by  the  sail  is 
very  localized  in  a narrow  wake  behind  the  sail  and  substantiates  the  assump- 
tion that  the  submarine  model  is  oriented  parallel  to  the  flow.  Figure  47 
for  X'/D  equal  to  10  shows  similar  results  except  that  the  profile  on  the 


side  of  the  body  has  started  to  depart  more  from  that  on  the  bottom 
opposite  the  sail. 

The  rake  was  modified  to  reduce  its  dimension  in  the  streamwise 
direction  to  a minimum  in  order  to  make  observations  as  close  to  the  pro- 
peller as  possible.  This  resulted  in  runs  at  X'/D  = 11.17  which  is  only 
about  2 inches  ahead  of  the  6 inch  diameter  propeller.  Some  results  are 
shown  in  Figure  48  where  the  top  to  bottom  asymmetry  produced  by  the 
wake  of  the  sail  is  still  clearly  evident. 

Figures  49  and  50  show  the  overall  effects  of  the  sail  on  the 
boundary  layer.  Figure  49  contains  all  the  plots  for  e equal  to  0 degrees 
which,  except  for  the  profile  at  X'/D  equal  to  1,  illustrate  the  effects 
of  the  sail . 

It  was  felt  that  the  flow  far  downstream  of  the  sail  might  separate 
were  it  not  for  possible  flow  entrainment  induced  by  the  propeller.  Any 
analysis  which  does  not  include  this  influence  on  the  flow  is  thus  in- 
complete. Figures  51  and  52  are  velocity  profiles  at  X'/D  equal  to  10  and 
11.17,  respectively,  with  and  without  the  propeller.  The  velocity  is 
left  dimensional  in  order  to  display  any  velocity  increase  from  flow 
entrainment,  if  it  exists.  At  X'/D  = 10  the  flow  is  identical  in  both 
cases  with  and  without  the  propeller,  indicating  no  flow  change  due  to 
entrainment,  but  it  is  apparent  that  the  propeller  is  starting  to  affect 
the  flow  at  11.17  diameters.  The  effect  is  more  pronounced  nearest  the 
body  where  the  velocity  difference  is  about  5 ft. /sec.,  an  increase  of 
9.0%,  whereas  in  the  freestream  the  increase  is  about  2 ft. /sec.,  a 1.7% 
increase  in  velocity. 
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V.  DISCUSSION 


A comprehensive,  experimental  program  of  study  considering  the 

flow-field  around  and  behind  a slender,  stern-propeller-driven  body  in 

5 

a wind  tunnel  at  nominal  ReQ  = 4.4  x 10  has  been  conducted.  The  work 
was  divided  into  four  tasks.  The  first  task  dealt  with  the  influence 
of  pitch  angle  variations  and  a sail-like  appendage  on  the  wake  behind 
a single  propeller  model.  The  second  task  dealt  with  the  wake  behind 
a dual  (side-by-side)  propeller  model  at  zero  pitch  angle.  The  third 
task  concerned  measurements  of  the  boundary  layer  over  the  body  itself 
and  the  influence  of  a propeller  on  the  body  layer  near  the  stern. 

The  fourth  task  involved  measurements  of  force  and  moment  coefficients 
for  both  models  over  a range  of  pitch  angle  from  -6°  to  +6°.  The  major 
results  and  conclusions  for  each  task  were  described  in  each  separate 


section. 


The  results  suggest  some  useful  areas  for  further  study.  First, 
the  influence  of  variations  in  propeller  parameters  such  as  the  number 
of  blades  and  solidity  should  be  determined.  It  will  also  be  interesting 
to  study  the  effects  of  roughness  for  boundary  layer  tripping  on  the 
blades.  Second,  since  the  sail  clearly  has  a large  influence  on  the 
flow,  a parametric  investigation  of  the  effects  of  this  and  other 
appendages  should  be  worthwhile.  Third,  further  studies  of  the  body 
boundary  layer  especially  nearer  to  the  propeller  are  sure  to  be  enlightening. 
Last,  the  whole  matter  of  making  good  turbulence  measurements  (including 
shear)  in  the  high  gradient  region  behind  a propeller  must  be  thoroughly 
considered. 
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Fig.  3(b)  SINGLE  PROPELLER  MODEL 


Fig.  4 CUTAWAY  OF  CYLINDRICAL  CENTER  BODY 
SHOWING  STRAIN  GAGE  BALANCE 


Fiqure  7:  Schematic  Of  the  Boundary  Layer  Rake 


Vertical  Traverse 


Horizontal  Traverse  Across 
Sail 


Horizontal  Traverse 


a)  Single  Propeller  Model 


Vertical  Traverse 


Off-Center  Vertical  Traverse 
] (Through  Hub) 


Horizontal  Traverse 


b)  Dual  Propeller  Model 
FIG.  Q WAKE  TRAVERSE  PATHS  (Looking  Upstream) 
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Fig.  16  Rolling  Moment  vs.  Pitch  Angle 

for  Both  Models  in  a Self  Propelled 
Condition 
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ANGLE  (DEG) 


Fig.  18  MEAN  AXIAL  VELOCITY  , FLOW  ANGULARITY,  AND  STATIC  PRESSURE 
HORIZONTAL  PROFILE  AT  Z/D=2,  Y/R=0,  ae=0° 
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Fig. 21  MEAN  AXIAL  VELOCITY,  FLOW  ANGULARITY,  AND  STATIC  PRESSURE 
HORIZONTAL  PROFILE  AT  Z/D=40,  Y/R=0,  ae=0° 


Fig.  22  MEAN  AXIAL  VELOCITY,  VERTICAL  PROFILES  AT  Z/D=40,  X/R=0 
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Fig.  33  MEAN  AXIAL  VELOCITY , FLOW  ANGULARITY,  AND  STATIC  PRESSURE 
VERTICAL  PROFILE  AT  Z/D=10 , X/R=0.  «e=-2° 
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IN  SAIL  WAKE  AT  Z/D=40,  Y/R=2.33 


MAXIMUM  AXIAL  TURBULENCE  INTENSITY  VS  Z/D 


FIG.  42  AXIAL  VELOCITY  DISTRIIUTIONS  FOR  THE  DUAL 
PROPELLER  MOOa  AT  Z 10  -10 


APPENDIX  A 


TABULATED  DATA  FOR  WAKE  SURVEYS 


WITH  SINGLE  PROPELLER  MODEL 
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B-3  (Continued) 
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TABLE  B-5  (Continued) 
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